Observation of Hybrid Soliton Vortex-Ring Structures in Bose-Einstein Condensates 
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We present the experimental discovery of compound structures comprising solitons and vortex 
rings in Bose-Einstein condensates (BECs). We examine both their creation via soliton- vortex 
collisions and their subsequent development, which is largely governed by the dynamics of interacting 
vortex rings. A theoretical model in three-dimensional (3D) cylindrical symmetry is also presented. 
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Quantized vortices have long been considered a signa- 
ture of superfluidity in Helium II P] . The recent emer- 
gence of Bose-condensed alkali gases provides a new and 
different arena for studying quantum fluids, their exci- 
tations, and interactions of the latter. In addition to 
observing vortices alone |^ Q , experimentalists have 
reported observations of grey solitons 5, 6] and their 
decay into vortex rings (VR) 0, ll| via the Kadomtsev- 
Petviashvili, or "snake" , instability 

In this Letter, we present the observation of a new class 
of topological excitations in BECs: hybrid structures 
originating from a combination of solitons and VRs. By 
'soliton', we mean a non-dispersive nonlinear wave whose 
phase profile and velocity are determined by its ampli- 
tude. Solitonic wavefronts extend on 2D sheets whereas 
vortex structures have a depleted core and phase singu- 
larity along a ID filament, which may loop to form a 
VR. We detail how these nonlinear excitations initially 
collide and how the subsequent compound structure dy- 
namics diverge from those of the superposed evolution 
of individual excitations. To date, the interplay between 
solitons and vortices has only been explored theoretically 

In previous work [Sj , we studied grey solitons and vor- 
tices generated from the quantum analog of shock waves 
by creating a single, narrow density defect in a BEC. We 
presently optimize the likelihood of collisions between 
these excitations by creating two symmetricall5M3laced 
defects via a mechanism similar to that used in [Sj. 

A thin wire blocks the 'coupling' beam from illuminat- 
ing the center of a magnetically trapped BEC (Fig. 1(a)), 
to create a symmetric double light-roadblock. We then 
inject two counter-propagating Gaussian-shaped 'probe' 
pulses along the symmetry axis of the trap. The pulses 
compress spatially b y se veral orders of magnitude via the 
slow-light technique Il4| , and ultimately fit entirely inside 
the BEC. The probe beams are resonant with the atoms' 
internal |1)-|3) transition while the coupling beam is res- 
onant with the |2)-|3) transition (Fig. 1(b)). 

When the pulses arrive at the wire's shadow, their 
propagation is effectively arrested and they are further 
compressed, as their speed depends on the coupling beam 
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FIG. 1: Defect creation, (a) Double light-roadblock set- 
up, (b) Energy level diagram; Mf refers to the projection 
of the atoms' angular momentum along the magnetic field 
quantization axis, z, of the 4-Dee electromagnetic trap K^ . 
(c) Density plots of |1) atoms from 2D simulations of defect 
creation using a 60 nm wire mask; initially, the BEC has 
4.0x10^' atoms. Light pulse intensities and durations are as 
in the experiment. 



intensity that rapidly falls to zero at the borders of this 
region ^J. Subsequently, atoms at the halted pulses' loca- 
tions are ejected, having first been transferred from initial 
state, |1), to |2) (untrapped) via a slow-light induced dark 
state. This results in two deep, narrow density defects in 
|1), symmetric about the BECs center (Fig. 1(c)). 

For the current experiment, we condense 3.6x10® 
sodium atoms in |1) and employ magnetic trap frequen- 
cies W2=27rx21 Hz and uj.j.—ojy—3.0LLJz- We work at tem- 
peratures well below the transition temperature to BEC 
so there is virtually no non-condensed component. The 
couphng (peak probe) Rabi frequency is 27rxl5.3 (3.8) 
MHz, and duration (probe 1/e half- width) is 12 (2.5) /iS. 
We use a 35 fim. wire mask. 

After creating the defects, we vary the duration of con- 
densate evolution in the magnetic trap, Attrap, and then 
release the condensate and let it expand for Atexp—^Q-Q 
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ms before using a slice technique [3, ll5| to image the 
central 20-30 /xm thick slab of the expanded cloud in the 
x-z plane. Figure 2(a) shows the 'control' experiment 
in which a single probe pulse is injected into the BEC, 
while Fig. 2(b) depicts the double-roadblock experiment 
in which interactions between excitations are evident. In 
both cases, the data are highly reproducible. 

The early dynamics that ensue from the creation of 
a single defect are detailed in Fig. 3 of 0. The defect 
splits into two density dips propagating at the sound 
speed in the +z and —z directions, while the back edge of 
each steepens due to the density dependence of the sound 
speed. In a classical fluid, this would lead to shock wave 
formation i.e. the back edges would develop infinite slope. 
In a superfluid, significant density changes over distances 
shorter than the healing length cannot occur . 

Consequently, as the density dips travel, they shed grey 
solitons each time their back edges become too steep [3 . 
As observed in Fig. 2(a), the solitons, although stable in 
ID, kink and decay into vortex structures via the snake 
instability, seeded predominantly by the transverse den- 
sity variation of the condensate depicted in Fig. 1(c). 

We now describe the experimental data in Fig. 2, com- 
paring and contrasting the double-defect images with the 
control experiment, before considering a corresponding 
theoretical simulation. Focusing mostly on structures be- 
tween the two initial defects, we see at Attrap=0.1 ms in 
Fig. 2(b) that solitons (pale stripes) have already shed 
from density dips originally emanating from the defects. 
High-density bands on the top and bottom result from 
the rush of fluid into the defect locations. The first soli- 
tons shed are the deepest, hence slowest moving, and 
therefore bound the regions of high density. In the cen- 
ter, between these first solitons, secondly shed ones have 
kinked so that they overlap on the sides but not in the 
middle, forming a checkerboard pattern. This is clearly 
seen at Attrap=0.4 ms where the images also show VR 
cores in the 'four corners'. Figure 2(b) at 0.1 and 0.4 
ms and Fig. 2(a) at 0.5 ms are comparable. The single- 
defect image shows additional shallower solitons in the 
wake of the travelling density dip, that are initially hid- 
den in the double-defect version by the mirror image of 
higher-contrast structures. At A<trap=l.l ms (Fig. 2(b)), 
soliton curvature has increased even more. Like at 1.3 
ms of Fig. 2(a), the deepest solitons form double-cusped 
shapes ('W on the top,'M' on the bottom) that span the 
cloud from side to side. The second pair of sohtons prop- 
agate toward one another further, closing off a region of 
fluid in the center. 

We have explored the 3D structure of the excitations 
by varying the y-position of the imaged slice to find that 
the system is highly cylindrically symmetric. Thus, the 
white loop in the center at 1.1 ms and the pairs of flat- 
tened loops seen in the shots corresponding to Attrap=l-3 
and 1.7 ms represent roughly ellipsoidal shells of low den- 
sity, a feature clearly absent in the single-defect case. 



Comparing the latter two time-points with Attrap=2-1 
ms in Fig. 2(a), we conclude that the centers of the deep- 
est solitons have decayed into VRs, forming the equato- 
rial components of their corresponding closed low density 
shells in Fig. 2(b). The upper and lower 'hemispheres' of 
the shells all derive from solitons formed after the first 
(deepest) ones that produced VRs. 

At Attrap=2.3 ms in Fig. 2(b), the closed structures 
appear to interlace and subsequently form a single low 
density shell at the center of the condensate (2.8-6.0 ms). 
This provides very strong evidence that the double-defect 
dynamics are not merely a superposition of mirrored 
single-defect dynamics. At Attrap=^-^ rns. Fig. 2(a) dis- 
plays a soliton fragment in the center with a VR below. 
At 4.4 ms in Fig. 2(b), the low density structure is closed 
and could not be formed from the superposition of ele- 
ments from the corresponding time point in Fig. 2(a). De- 
spite being larger at Attrap=5 and 6 ms (Fig. 2(b)) than 
in previous frames, this compound structure is temporar- 
ily stable. At 6.6 ms, we observe that what was a closed 
shell has sheared and come apart at the sides. The re- 
sulting structures break cylindrical symmetry, highlight- 
ing their sensitivity to torque from the laser beams, ev- 
ident in Fig. 1(c). Furthering the case for interactions 
rather than superposition, the last two frames of Fig. 2(a) 
have a long-lived VR in the center, in stark contrast with 
Fig. 2(b) at 6.6 ms, which has segments that could not 
be constructed from VRs alone. 

We measure the soliton and shell half- widths to range 
from 3 to 6um, comparable to a calculated healing length 
of 4 /im [l2|. We also note the collective modes in the 
BEC produced by the initial rush of fluid into the defects' 
locations. The high density bands in both Figs. 2(a) and 
2(b) bend, contract, and expand. We find these modes 
are highly dependent on the initial defect shape. 

To complement our observations, we have performed 
simulations based upon a generalized Gross-Pitaevskii 
(CP) description ^ for three phases of BEC evolution 
as follows: We first simulate defect creation via slow- 
light propagation on a 2D spatial grid, as in Fig. 1(c). 
The atoms in |2) are quickly ejected from the trap and 
we map a symmetrized version of the resulting defect in 
the |1) condensate onto a 3D cylindrically symmetric grid 
[l9| . We then calculate in-trap evolution of the |1} atoms 
for a duration Attrap and finally continue to follow the 
BEC evolution after the trapping potential is turned off. 

The results of simulations for Atexp=^'^ ms with vary- 
ing Attrap are shown in Fig. 3. (After 14 ms of cloud 
expansion, all topological dynamics are frozen out and 
further evolution simply leads to spatial magnification 
of the developed features; the BEC's aspect ratio differs 
from that at 19.9 ms by ~15%.) We see the same promi- 
nent features in the simulation as in the experiment, al- 
lowing us to identify the nature of the low density shells. 
Phase information confirms the existence of solitons and 
identifies the creation and annihilation of VRs. As we 
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FIG. 2: Experimentally observed resonant transmission images of BECs (transmitted laser intensity normalized to incident) 
after they are illuminated with coupling and probe pulse(s). In (a) ((b)), one (two) probe pulse(s) is (are) hitting the cloud, 
there are initially 4.0 (3.6) x 10® condensed atoms, and the imaged slice is 30 (23) ^m thick. Times cited refer to evolution time 
in trap, Aftrap, after defect creation; Atea:p=19.9 ms. Arrows indicate initial defect planes. 
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FIG. 3: Simulation for Ate2:p=14 ms (Aft,ap on plots) with 
parameters the same as in Fig. 1(c). To compare with experi- 
ment, we plot e"^"'^'^', where Dy{x, z) is the optical density 
along y over the central 22 //m slab of the BEG. 

detail below, the simulations indicate that the expansion 
'restabilizes' solitons that would have been subjected to 
the snake instability, if kept in the trap. Experimental 
parameters have been chosen to cover a regime in which 
VRs are produced but parent solitons are not destroyed. 

The simulation follows that of the experiment up to 
ms, after which point evolution is slightly ac- 
celerated. Other discrepancies include reduced soliton 
curvature and less dramatic collective mode breathing 
than in the experiment. This could be a result of the 
different initial defect separation used in the theory to 
generate a sequence most comparable to the data. The 
simulation's cylindrical symmetry is an added approxi- 
mation, but it generally describes the observations well. 

Despite these differences, the simulations provide us 



with the history of the experimentally observed struc- 
tures that must be expanded to be optically resolved. 
Calculations show that if two solitons are incident on 
one another, the time that it takes them to reach each 
other and their degree of kinking at that point have dif- 
ferent density dependencies. During expansion but be- 
fore interaction energy is completely converted to kinetic, 
soliton translational velocity and the rate of soliton kink- 
ing both decrease as BEC density drops, but at different 
rates. These different density dependencies do not only 
imply that expansion evolution is highly important and 
different than in-trap evolution ; they also lead to the 
notion that expansion sequences succeeding two different 
in-trap durations are not mere duplicates offset from each 
other in time. We see variations in such sequences, typ- 
ically manifested by the occurrence of VR creation and 
annihilation at slightly different points in the evolution. 

As an example, we present the simulated release dy- 
namics for Aftrap=3.0 ms in Fig. 4. As in the experiment, 
the first soliton (blue) shed from each of the waves travel- 
ling toward the EEC's center is the deepest, slowest, and 
least stable, leading to decay into VRs at Atea;p=l-5 ms. 
Solitons created after the first ones (green through red, in 
order of creation from top defect) propagate across the 
center line of the BEC (Aiea;p=1.5 and 3.0 ms); those 
with opposite propagation directions sometimes overlap. 
As the third and fourth generated solitons reach the 
newly formed VRs on the side opposite whence they orig- 
inated, we observe pairs of low density shells (Aiea:p=4.0 
ms). By examining the phase, we confirm, as discussed 
for images 1.3 and 1.7 ms of Fig. 2(b), that these shells 
are made from VRs and transiently passing soliton fronts. 

The second deepest (secondly shed, green) solitons 
cross in the center of the BEC as their midpoints bend 
toward one another (also Aiea;p=4.0 ms). As these fronts 
reach the VRs at 5.0 ms (again forming a pair of low 



4 
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FIG. 4: Simulation of soliton VR dynamics during expansion after Attrap~3.0 ms. Times cited are Atexp- Only solitons 
propagating in —z direction (downward) from the +z (top) defect are shaded (blue through red, in order of creation), to assist 
in tracking the movement of each from frame to frame. Dots (crosses) indicate vortex cores for z>0 only with positive (negative) 
circulation with respect to the plane of the page. The vortex colors correspond to parent soliton colors, although the green 
cores are generated by the mirror image of the correspondingly shaded soliton. Density ranges from to {7.5, 4.5, 4.1, 3.6, 1.3, 
0.92, 0.42} X 10^^ cm~"^ for each of the times presented, respectively. Simulation parameters are the same as in Fig. 1(c). 



density shells), the central portions of the sohtons curve 
further, producing within themselves two smaller VRs 
with opposite circulation from the larger ones. These 
truly hybrid low density structures, consisting of par- 
ent soliton segments and VRs, 'dislocate' from the rest 
of the soliton fronts. They each comprise two concen- 
tric VRs embedded in a roughly hemispherical solitonic 
shell. The smaller VRs are evident from the wavefunc- 
tion's phase even though they are not distinguishable in 
the density plots. The general motion of the hybrid seg- 
ments is largely determined by the velocity fields from 
the vortices present. Like an umbrella turning inside-out, 
the curvature of these hybrid structures reverses through 
nis, as the inner rings propel themselves away 
from one another and the outer rings propel toward the 
center line. This occurs in the experiment for Attrap='2-3 
ms and in Fig. 3 at 1.7 ms. The VRs behave as the 'phan- 
tom' propellers in the compound structures. 

As the outer VRs propagate toward each other, 
a closed shell forms from the two hybrid segments 
{Atexp—Q-5 ms). It is an ellipsoidal solitonic shell incor- 
porating four VRs of alternating circulation that encircle 
the structure like lines of constant latitude on a globe. 
The outer VRs eventually annihilate (quantized circula- 
tion at the VR cores disappears). Subsequently, the two 
segments that formed the central shell lose their curva- 
ture, and the shell breaks up along 2=0. This is seen at 
Aiea:p=14 ms (a close-up of the BEC at Attrap—3-0 ms 
in Fig. 3) and should be compared to the experimentally 
observed shearing at 6.6 ms in Fig. 2(b). 

In light of these observations it is natural to ask 
whether low density shells can exist as stationary states 
in a BEC. We have shown that the GP equation ex- 
hibits stationary solutions in spherical symmetry consist- 
ing of an inner sphere of almost constant density sur- 
rounded by low density shells close to Bessel-function 
type spherical standing waves sin(fcr)/r known from the 
linear Schrodinger equation. This has been indepen- 
dently treated by others The infinite space solu- 
tions have an infinite number of particles and concentric 
shells. In a trap, similar solutions with finite particle 



number and number of shells are possible; details will 
be given elsewhere. Unlike the observed structures dis- 
cussed above, the spherical stationary states show no vor- 
tex core structures. However, the observed hybrid shells 
might play a role in the decay of the stationary solutions. 

In brief, we have presented the experimental observa- 
tion and theoretical confirmation of low density shells 
in a fully 3D BEC, consisting of complex hybrid soliton 
vortex-ring structures. VR propulsion, attraction, and 
annihilation heavily influence the structures' dynamics. 
The BEC's decreasing density, after a confining trap is 
switched off, contributes critically to the dynamics, as 
the resulting regime is on the cusp of stability of solitons 
typically subject to the snake instability in 3D. 
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